Simulation of the spin-boson model with superconducting phase qubit coupled to a transmission line 
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Based on the rapid experimental developments of circuit QED, we propose a feasible scheme to simulate a 
spin-boson model with the superconducting circuits, which can be used to detect quantum Kosterlitz-Thouless 
(KT) phase transition. We design the spin-boson model by using a superconducting phase qubit coupled with a 
semi-infinite transmission line, which is regarded as bosonic reservoir with a continuum spectrum. By tuning the 
bias current or the coupling capacitance, the quantum KT transition can be directly detected through tomography 
measurement on the states of the phase qubit. We also estimate the experimental parameters using numerical 
renormalization group method. 
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Quantum simulation is one of the original inspirations for 
quantum computing, proposed by Feynman |10 to solve the 
difficulties of simulating quantum systems on a classical com- 
puter. Simulating an arbitrary quantum system by the most 
powerful classical computer is very hard for large scale quan- 
tum systems, because of the exponential scaling of the Hilbert 
space with the size of the quantum system. Quantum phase 
transitions (QPT) ^ at zero temperature play a key role in 
the occurrence of important collective phenomena in quan- 
tum many-body systems, which occur as a result of compet- 
ing ground state phases. Similarly, simulation of QPT with 
classical computer is also difficulty since it is usually relevant 
to a quantum many-body system. 

On the other hand, the spin-boson model a two-level 
system linearly coupled to a collective of harmonic oscil- 
lators, is a typical model to study decoherence effects and 
QPT. Those dissipative spin systems 14]] are very interest- 
ing because they display both a localized (classical) and de- 
localized (quantum) phase for the spin. Spin-boson model 
has been primarily investigated by numerical renormalization 
group (NRG) m la] and find that different spectral functions 
of the bath may induce various kind of QPT. Many efforts 
have been made to observe such environment-induced QPT in 
various systems, such as mesoscopic metal ring |71], single- 
election transistors lES] and cold atoms ifioll . etc.. 

In this paper, we propose a feasible scheme to simulate the 
spin-boson model with superconducting circuits, which can 
be used to observe the notable quantum KT phase transition. 
Here we focus on the ohmic case (s = 1), which can be 
mapped on the anisotropic Kondo model. The model shows 
a Kosterlitz-Thouless (KT) quantum phase transition Jsj], sep- 
arating the localized phase at a > Ofc from the delocalized 
phase 2X a < Uc lUHl , where a represents the strength of 
the dissipation and ac is the critical value. Our idea is in- 
spired by previous efforts to study the superconducting cir- 
cuits as an artificial atom lfl2ijl8ll . In the paper we design the 
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no. 1 : A phase qubit capacitively coupled to a semi-infinite trans- 
mission line simulates the spin-boson model. 



spin-boson model using a superconducting phase qubit cou- 
pled with a semi-infinite transmission line, which is regarded 
as the qubit' s environment with gapless spectra. In this setup, 
the spectral function of the bosonic bath is treated as Ohmic. 
By tuning the bias current or the coupling capacitance, the 
coupling between the spin and the environment can be con- 
trolled. So the states of the qubit may transit from delocal- 
ized phase to localized state, which can be directly observed 
by measuring the phase qubit. We also estimate the experi- 
mental parameters of this transition with NRG method. Com- 
paring with other candidates of spin-boson modellfTUloll, the 
proposed experimental setup based on superconducting sys- 
tem lll2t - il4ll may have some distinct advantages, such as the 
parameters in the model are tunable through experimentally 
controllable bias cuiTents or driving microwaves, the bosonic 
reservoir with a continuum spectrum can be simulated easily, 
and the measurements are of high-fidelity, etc.. 

Our designed experimental setup is illustrated in FiglT] 
an approximately ID transmission line (with the length as 
L oo) is capacitively coupled to a current-biased phase 
qubit. The phase qubit is an artificial spin, while for rela- 
tively low frequencies the transmission line is well described 
by an infinite series of inductors with each node capacitively 
connected to ground and then can be considered as a boson 
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bath lllSll . The Hamiltonian of the whole system can be writ- 
ten as 



H = H„ + Hr, 



Hi, 



where 



Ej cos ( 



' 2 '27r""' 27r 
is the Hamiltonian of the current-biased phase qubit. 



is the Hamiltonian of the transmission line, and 
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(4) 



is the interaction between the phase qubit and the transmission 
line. In these equations Ej ~ is the magnitude of max- 
imum Josephson coupling energy, Iq is the critical current of 



the junction, $0 = 57 is the superconducting flux quantum, 
Ih is the bias current, (j>j is the phase difference of the junc- 
tion, Cj and Co are the junction and coupling capacitance, c 
and I are the capacitance and inductance per unit length, re- 
spectively; '(^(a;, = (7(2;', t)(ia;' is the collective charge 
variable on the transmission line. 

For the current-biased Josephson junction |19J, the charge 
operator Q = C^cjyj and phase difference operator (j)j have 

the commutation relationship [(/) J, Q] = 2ei. Quantum me- 
chanical behavior can be observed for large area junctions 
in which Ej ^ Ec = e^/2C when the bias current is 
slightly smaller than the critical current < /q. In this 
regime, the last two terms of the Hamiltonian of the phase 
qubit (as Eq.(|2|) can be accurately approximated by a cubic 
potential U{(j)j) parameterized by a barrier height AU (lb) = 
(2\/2/o(l>o/37r)[l - (/{,//o)]^/^ and a quadratic curvature at 
the bottom of the well that gives a classical oscillation fre- 
quency cjp{h) = 2i/4(2^/o/$oC)i/2[l _ (h/Io)]'/^ wifli 
the capacitance C = Cj + Cq. The lowest two of the quan- 
tized energy levels in the cubic potential are considered as the 
qubit states, {|0),|1)}, where uiq ~ 0.95cjp is the energy 
difference between the ground state and the first excited state 
of the phase qubit. The states can be fully manipulated with 
low- and microwave frequency control currents, which can be 
chosen as /t = Idc 
expressed as 



I^iw{t) fil9i1 . so the phase qubit can be 



nn — -Ox — cr2, 



(5) 



/^^ (i) and A 



where cr^.z are Pauli matrices, e = 
huj\Q. In the system, the charge operator can be described as 

The corresponding Euler-Lagrange equation for the trans- 
mission line is a wave equation 

\d'^^{x,t) d{x,t)S{x) I d^'d{x,t) 
c fe2 Co 2 dt^ 







with the mode speed v = I /Vic. This is analogous to the 
problem of the Schrodinger equation with a delta function po- 
tential. By separation of variables, the spatial part of the solu- 
tion of the modes is of the form 



\ A cos(kx) symmetric modes, 

17 g ix) — \ 

I A sin (/ex) antisymmetric modes, 
where Q < x < L. According to the boundary condition 

1 ai9(0) d{L) 



(6) 



c dx 



Co 



in the limit Co ^ c, the spatial mode can be expressed as 
k « 2^ , where m is odd for symmetric modes and is even 
for antisymmetric modes. The eigenfrequencies of the modes 
are w = kv. The time dependent part for transmission line 
still has the form of Euler-Lagrange equation 

4>m{t) + W„(/),„(i) = 0. 

From the above equation we can obtain the time dependent 
Hamiltonian for multi-mode transmission line 



(7) 



as a function of (j)m{t) and its canonically conjugate momen- 
tum p„ = /(/>™(t) with [0m, Pm'] = ihSrnni'- 

To diagonalize the Hamilitonian HT{t), we use the usual 
relations 



hbJmC 2L 



(8a) 



-tM - alit)] (8b) 



by introducing the bosonic creation and annihilation op- 
erators a: [armO-l^i] — 6mm'- Then the collective charge 
variable in the semi-infinite transmission line is ■d{x,t) = 

cos ^x (m odd), 



ni—1 



^-[amit) + alit)] 



2L 



sm^^x (meven). 
The voltage at the first end of the transmission fine (x — 0) is 



y(0,i) 



c dx 



E 

n=l 



[a^{t) + ai{t)l (9) 



where n = m/2 = 1,2, ...ric- Here m is even for the maximal 
voltage. The voltage is zero for odd m and thus this part can be 
neglected. Substituting Eq.(|9]) into Hint in Eq.(|4|i, we obtain 
the interaction Hamiltonian between the phase qubit and the 
transmission line given by 
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(10) 
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After rotating the frame of the phase qubit (cr^ a^TCTy ^ 
az), the Hamiltonian of the circuit without microwave fre- 
quency control current (e ~ 0) can be written as the standard 
spin-boson model, 

H = -^(^x + fiLLi„alan + Y X!'^"^"" + ^^^^ 



where 



Co 



2AhLU„ 
CLc 



are the coupling strengths between 



the spin and the modes of transmission line, and the frequen- 
cies of modes a;„ = It is notable that the coupling 
strengths can be controlled by the bias current (or the cou- 
pling capacitance Cq if it can be modified in the experiments), 
so the parameters in the spin-boson model realized in this su- 
perconducting circuits are tunable and thus the system is suit- 
able to be used to observe rich phenomena in the spin-boson 
model. 

The spectra is gapless when the length of the transmission 
line is sufficient large. In this case, the bosonic bath can be 
characterized by its spectral function isfc 



(12) 



In the proposed experimental setup, the obtained spectral 
function is actually Ohmic case s ~ 1. In additional, the di- 



c I — 
hbjy \/l/c mea- 



mensionless parameter a = ^ \/W^ - 
sures the strength of the dissipation, which is determined by 
the coupling strengths A„ and may be modified with the bias 
current and the coupling capacitance Cq. 

We now turn to address a method to observe one of the most 
interesting phenomena in the spin-boson model: the KT phase 
transition. In the dissipation case, the spin-boson model has a 
delocalized and a localized zero temperature phases separated 
by the KT transition at the critical value ac — 1 (for the 
unbiased case of e 0). In the delocalized phase at small 
dissipation strength a, the ground state is nondegenerate and 
represents a damped tunneling particle. For large a, the dissi- 
pation leads to a localization of the particle in one of the two 
az eigenstates, thus the ground state is doubly degenerate. 

In the proposed experimental setup, the strength of the 
dissipation is proportional to the controllable parameter as 
a oc A/27r, which can be changed by tuning the bias cur- 
rent /f,. The parameter a can be experimentally modified 
in the regime 0.2 ~ 3.0 if one takes the following typical 
data: the junction capacitance Cj = 0.85 pF, Co = 5Cj, 
the impendence of transmission line z = \/TJc = 50 fl, and 
A ~ fujjp/2. Therefore, the localized and delocalized states 
as well as the KT phase transition may be observed by tuning 
the parameter A/ 2 through the biased current I^. 

The different phases and the KT phase transition can 
be demonstrated through measuring the populations of the 
qubit states (ctz), which can be detected by a three- 
Josephson-junction superconducting quantum interference 
device (SQUID). The result of measurement can be defined 
as 5P = h\Pe ^ PgV where Pe{g) denotes the population of 
excited (ground) state of phase qubit. (5P = means that the 
system stands at the delocalized phase, while 5P = 0.5 stands 
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FIG. 2: The differences of the occupation probability SP — \Pe — 
Pg\/2 for several values of e. The sudden change of the populations 
characterizes the KT type quantum phase transition in the spin-boson 
model, which occurs from the delocalized phase at 5P = to the 
localized phase at SP = | . 



at the localized phase. The 6P as a function of the strength a, 
calculated with numerical renormalization group method, is 
plotted in FiglJ] The readout technique of phase qubit 1I20I1 
is achieved by applying a short bias current pulse SI{t) (less 
than 5 ns) that adiabatically reduces the well depth AU /hujp, 
so that the first excited state lies very near the top of the well 
when the current pulse is at its maximum value. It has been 
shown that only the expectation of cr^ in the phase qubit can be 
measured; however, the other direction measurements can be 
achieved by applying a transformation that maps the detected 
eigenvector onto the eigenvector of cr^ before detecting. The 
ratio of the tunneling rates for the two states |1) and |0) is 
about 200, so the fidelity of measurement of the phase qubit 
is about 96% when properly biased. 

Let us briefly introduce the critical behaviors of the 
model studied by NRG, which is an efficient way to treat 
the spin-boson model with a broad and continuous spec- 
trum of energiesiH 0]. The NRG method starts with a 
logarithmic discretization of the bosonic bath in intervals 
[A-("+i'a;c, A-"wc](?i = 0, 1, 2, . . .), where A > 1 is called 
as the NRG discretization parameter We take the surface 
plasma frequency as the cutoff energy {ujc = 10^^ Hz) and 
energy unit. After a sequence of transformations, the dis- 
cretized model is mapped onto a semi-infinite chain with the 
spin representing the first site of the chain. The spin-boson 
model in the semi-infinite chain form is diagonalized it- 
eratively, starting from the spin site and successively adding 
degrees of freedom to the chain. The exponentially growing 
Hilbert space in the iterative process is truncated by keep- 
ing a certain fraction of the lowest-lying many particle states. 
Due to the logarithmic discretization, the hopping parameters 
between neighboring sites fall off exponentially, going along 
the chain corresponds to accessing decreasing energy scales 
in the calculation. From the energy of the first excited state, 
the properties of the total system can be presented according 
to the parameters, i.e.. A, a, e. It should be mentioned that, 
when the controlling parameters are close to the critical val- 
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FIG. 3: Scaling of the flow of the many-particle levels i5jv(j') for 
fixed s = 1, e = 0, and A tuned in the regime 2.5 x 10^'^ ~ 3.75 x 
10""", and the corresponding a is changed in 0.5 0.75. The NRG 
parameters are Ns — 100, Nt = 6, and A — 2.0. 

ues, the numerical computation would occur small deviation 
by NRG. The result of the critical value that obtained directly 
by numerical calculation is less than practical one, which is 
analogous to that e is finite small value. In order to illuminate 
this errors, we consider the cases of finite small values e to 
show the effect on a, as described in Fig ID 

We expect to observe scaling behavior in all physical prop- 
erties for q;(A) — > ac(Ac). An example is shown in Figl3]for 
various values of q;(A). In this way we can easily determine 
the critical value Uc from the relation T* ~ const, x A^'^ cx 
/\i/iaa-a) ^jy* jg jj^g yalue of N wherc the first excited state 



reaches the value = 0.3, the NRG discretization pa- 

rameter is A = 2.0 in the paper). From the sets of data { A^, a} 
that approaching critical point in Figl3] the critical value can 
be determined as ac = 1.093 via extrapolation, which is con- 
sistent with the exact value etc = 1 in the small A limit 
The critical value separate two different phases, which can be 
observed directly by measuring the populations of the qubit 
states. 

Before ending the paper, we make two additional remarks. 
If the coupling capacitance Cq can be tuned, we can also de- 
tect the quantum KT phase transition through modifying Cq 
but fixing A. On the other hand, we just discussed the Ohmic 
spectra in the paper since the resistances of the transmission 
line are neglected; however, the sub-Ohmic spin-boson model 
^ can be designed if a RC-dominate transmission line is 
considered. In this case, the critical properties of the sub- 
Ohmic spin-boson model can be observed in a slightly ex- 
panded model. 

In summary, we have presented a scheme to reaUze the 
spin-boson model by using a phase qubit coupled to a semi- 
infinite transmission line. By tuning the bias current or the 
coupling capacitance, the notable KT phase transition from 
the delocalized phase to the localized phase can be directly 
observed through measuring the states of the phase qubit. We 
have also estimated the required experimental parameters by 
using numerical renormalization group. 
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